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APPLICATION FOR PATENT 

Inventors: Hedva Spitzer and Eilon Sherman 

Title: A method for automatic color and intensity contrast adjustment of still 
and video images 

5 FIELD AND BACKGROUND OF THE INVENTION 

The present invention relates to image processing and, more particularly, to an 
automatic method for correcting the contrast and intensity of color still and video 
images. 

The perceived colors of visible objects are determined, both by the physics of 
1 0 light reflection and by the way in which the visual system processes the reflected light 
that it receives. With regard to the physics, the physical color of the light reflected to 
the visual system by visible objects is determined, both by the reflectivity spectra of 
their surfaces, and by the color spectrum of the illuminating light. Photographs of 
scenes, including both still pictures and motion (e.g. video) pictures, whether recorded 
15 by analog means (photographic film) or digital means (video cameras), normally are 
perceived differently from the way the scenes themselves would be perceived by direct 
vision. 

Most prior art methods dealing with color image processing of still and video 
pictures, and specifically with color and intensity contrast, approach the subject of 
20 improving or adjusting such contrast and/or intensity using principles unrelated to the 
human visual system and its physiology. A previous patent awarded to one of the 
present co-inventors, U.S. Patent No. 5,771,312, describes physiologically-based 
color constancy adjustment of pictures, but does not deal with either color or intensity 




contrast. We are not aware of any method similar to the one proposed herein, that 
aims to emulate the human vision for improving color contrast and intensity. Present 
color contrast adjustment methods thus lack many of the physiological-like features 
that can lead to much improved color pictures. 

5 There is thus a widely recognized need for, and it would be highly 

advantageous to have, a satisfactory method for performing color contrast correction, 
either automatically or interactively, on both still and video (motion) pictures. There 
is also a widely recognized need for, and it would be highly advantageous to have, a 
satisfactory method for performing local intensity correction, either automatically or 
10 interactively, on both still and video pictures. 

SUMMARY OF THE INVENTION 

The underlying concept of the present invention is to process color pictures in a 
manner similar to that in which the neurons of the visual system process signals 
related to color vision to achieve color contrast, and to adjust intensity, which is a 

1 5 component of color. 

According to the present invention, there is provided a method for correcting the 
color contrast of a scene, the scene including an intensity spectrum at each of a 
plurality of pixels, the method comprising the steps of: a) providing a red image, a 
green image, and a blue image, each image having a pixel value at each of the 

20 plurality of pixels; at each pixel: b) computing a center red response, a center green 
response and a center blue response based on the images; c) computing a surround red 
response, a surround green response and a surround yellow response based on the 
images; d) computing a red, a green and a blue on-center non-filtered and filtered 
opponent response, based on the center and surround responses; e) computing a red, a 
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green and a yellow off-center non-filtered and filtered opponent response based on 
the center and surround responses; f) computing a red, a green and a blue double- 
opponent (or "do") response and a corresponding filtered do-response based on the 
on-center and off-center non-filtered and filtered opponent responses; g) computing a 
5 red, a green and a blue do-remote signal based on a set of responses selected from the 
group consisting of the on-center filtered opponent responses and the filtered double- 
opponent responses; and h) for each pixel: correcting each of the red, green, and blue 
double-opponent responses for color contrast using respectively the red, green and 
blue do-remote signals, thereby producing new red, green and blue double-opponent 
10 responses. 

Preferably, the method for correcting the color contrast of an image further 
comprises: at each pixel: i) computing a yellow non-filtered and filtered center 
response and a non-filtered and filtered blue surround response; j) computing a 
yellow non-filtered and filtered double-opponent response based on the non-filtered 

15 and filtered yellow center and blue surround responses; k) computing a yellow do- 
remote signal based on a response selected from the group consisting of the filtered 
yellow center response and the filtered yellow double-opponent response; and 1) for 
each pixel: correcting the yellow double-opponent response for color contrast using 
the yellow do-remote signal, thereby producing a corrected yellow double-opponent 

20 response. 

Preferably, the method for correcting the color contrast of an image further 
comprises: for each still image and for each pixel: a) inversely transforming the red, 
green and blue corrected double-opponent responses to obtain corresponding new red, 
green and blue opponent responses; and b) transforming the new red, green and blue 
25 opponent responses into corresponding new red, green and blue center responses. 



3 




Preferably, the method for correcting the color contrast of an image further 
comprises: for each picture and for each pixel: a) inversely transforming the corrected 
yellow double-opponent response to obtain a new yellow opponent response. 

The method of color contrast adjustment according to the present invention can 
5 also be used to decrease artifacts of image compression (especially JPEG pictures), 
and to decrease artifacts of both compressed and uncompressed picture transmission 

According to the present invention there is further provided a method for 
adjusting an achromatic contrast of a scene, the scene including an intensity spectrum 
at each of plurality of pixels, the method comprising the steps of: a) providing an 
10 image that has an intensity value at each of the plurality of pixels; at each pixel: b) 
obtaining an adapted opponent center response using a plurality of the pixel intensity 
values, and c) at each pixel, correcting the achromatic contrast using the adapted 
opponent center response. 

Preferably, the method for adjusting an achromatic contrast of a scene further 
15 comprises obtaining an adapted opponent surround response, wherein the step of 
correcting for achromatic intensity contrast includes subtracting the adapted opponent 
surround response from the adapted opponent center response. 

According to the present invention there is further provided yet another 
method for adjusting an achromatic contrast of a scene, the scene including an 
20 intensity spectrum at each of a plurality of pixels, the method comprising the steps of: 
a) providing an image having a pixel value at each of the plurality of pixels; at each 
pixel: b) computing a center response based on the image; c) computing a surround 
response based on the image; d) computing a non-filtered and filtered on-center 
opponent response, based on the center and surround responses; e) computing a non- 
25 filtered and filtered off-center opponent response based on the center and surround 
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responses; f) computing a double-opponent response and corresponding filtered 
double-opponent response based on the filtered on-center and off-center opponent 
responses; g) computing a do-remote signal; and h) for each pixel: correcting the 
double-opponent responses for achromatic contrast using the do-remote signal. 
5 The method of achromatic contrast adjustment according to the present 

invention can be used to decrease artifacts of image compression (especially JPEG), 
to enhance achromatic contrast in pictures taken under poor lighting conditions or 
with poor equipment, and to decrease artifacts of both compressed and uncompressed 
image transmissions. Both color contrast and achromatic contrast adjustments as 

10 described herein can be applied to any type of arrayed data (or image), for example to 
an image represented in polar coordinates, etc. 

Implementation of the algorithm (methods) for both color and intensity 
contrast correction, and also for object recognition, of the present invention may take 
one or more of this preferred forms: 

15 1. a) Hardware: i. implementation on a dedicated chip; ii. implementation on a 
dedicated computer slot; iii. integration with another chip; and iv. integration with any 
other chip computer slot. 

b) Software: i. implementation as an add-on software to the hardware devices 
of 1(a); ii. implementation as a stand-alone software package; and iii. integration with 

20 another software package (such as CAD, image processing, etc.). 

2. Enhancement of color contrast and intensity for still and video cameras. 
Cameras may include digital cameras and miniature video cameras (for video over IP 
applications). 

3. Enhancement of color contrast and intensity in any imaging application, such 
25 as: a) Biomedical - endoscopy, diagnostic imaging (such as CT, MRI, etc), b) Science 
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and research - any imaging for astronomy, chemistry, physics, biology, computer 
science, etc. c) Industrial - print industry (including digital-print), inspection, artificial 
imaging and robotics of any kind, automatic character or object recognition (OCR or 
ATR). d) Military and security applications - satellite images, long and short distance 
5 imaging and surveillance (ground based and aerial), displays, night vision and thermal 
vision (as false colors), e) Sports, news, fashion, commercials and advertisements, f) 
Professional, amateur and artistic photography. 

4. Communication: a) Elimination of artifacts in transmitted images (still and 
video) due to compression (such as JPEG and MPEG) and transmission media affects. 

10 Applications may include video conferencing, broadcast, TV (of any format), HDTV. 

5. Software applications: a) Enhancement of color contrast and intensity of 
images in a privately owned web-site, b) Same as (a), but as a download software, c) 
Same as (a) as a stand-alone software package, d) Same as (a) as a service (or even as 
a link) to another web site. End users may come from one or more of the the fields of 

1 5 interest of 3 and 4. 

The underlying concept of the present invention is to process color pictures in a 
manner similar to that in which the neurons of the visual system process signals related 
to color vision to achieve color contrast. Ideally, the input to the present invention is the 
intensity spectrum of a scene as a function of wavelength, measured at each pixel in an 

20 array of pixels. The term "scene", as used herein, refers to a set of objects that produce, 
by emission, by scattering or by reflection, visible light that can be perceived by a 
viewer. For example, a person viewing a pastoral scene sees light reflected by objects in 
the scene such as trees, grass, rocks and soil, as well as the scattered blue light of the sky. 
In each of the viewer's eyes, light from the scene is focused by the lens of the eye onto 

25 the retina of the eye. The light impinging on each point of the retina is characterized by 
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an intensity spectrum that spans the full range of wavelengths of visible light. In the case 
of an imaging apparatus such as a digital camera, the optics of the imaging apparatus 
typically focuses the light onto an array of sensors such as charge coupled devices; each 
sensor then produces a set of signals that correspond to spectral subbands of the light 
5 from one pixel of the scene. Ideally, the intensity spectrum would be multiplied by the 
spectral response function of each of the types of photoreceptor cells of the retina (red 
cones, green cones, and blue cones) to incident light, and integrated with respect to 
wavelength, thereby providing, at each pixel, a red intensity value, a green intensity 
value, and a blue intensity value. Collectively, the red values, the green values, and the 
10 blue values are examples of what is referred to herein as "images": rectangular arrays of 
values, one value per pixel. These values then are processed according to the algorithm 
of the present invention to provide images corrected for color contrast and intensity 
contrast. 

This ideal input rarely is attainable. Therefore, the scope of the present invention 
1 5 includes the processing of images obtained by other means. For example, the three input 
images may be in the form of analog signals from transducers whose spectral responses 
are similar to the spectral responses of cone cells, in which case the intensity values are 
electrical signals, typically voltage levels. These analog signals may be processed 
directly, using an embodiment of the present invention in analog hardware. 
20 Alternatively, these analog signals may be digitized, and processed digitally according to 
the present invention. Usually, however, the input to the present invention consists of 
digital images, such as are acquired by video cameras, that come ultimately from 
transducers whose spectral responses does not necessarily match the responses of cone 
cells. In that case, the digital pixel intensity values can be transformed to photoreceptor 
25 response coordinates, or "fundamentals", corresponding to the spectral responses of the 
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three types of cones, the responses of similar transducers that are used as "cones", or 
inputs from any other color scale.. 

The most common color coordinate system for digital color images is the so- 
called red-green-blue, RGB, or chromaticity, coordinates. Digital images in other three- 
5 color schemes, such as yellow-cyan-magenta, may be transformed mathematically to r- 
g-b. The transformation from RGB coordinates, or from CIE Yxy coordinates, to 
photoreceptor coordinates may be found, for example, in G. Wyszecki and W. S. Styles, 
"Color Science" (Wiley, 1982), pages 139 and 615. In what follows, all references to 
"red", "green" and "blue" will be to photoreceptor response coordinates, and not to 
10 chromaticity coordinates, and instead of the RGB scale, use will be made of the LMS 
(long-medium-short wavelength) scale. When "red" "green" and "blue" are specified, 
they are used as equivalents to L, M, and S respectively. 

The present invention includes an algorithm in the spirit of that presented by Ronen 
Dahari and Hedva Spitzer in an article titled "Spatiotemporal adaptation model for 
15 retinal ganglion cells", published in the Journal of the Optical Society of America Series 
A, Volume 13 Number 3 (March 1996), which article is incorporated by reference for all 
purposes as if fully set forth herein. The paper by Dahari and Spitzer presents a model 
for the adaptation of visual perception to changing intensity of illumination. It has been 
conjectured that color contrast perception and intensity contrast perception work by an 
20 analogous mechanism, with modifications as described herein. 

V4^*^ / Flt3w 1 is a schematic description of a cross section of the human retina, 
snowii 



lowing that the retk^a consists of five layers of cells, primary receptors 4, horizontal 
cells 6, bipolar cells 8, amafciijie cells 10, and retinal ganglion cells 12. Receptors 4 
shown here are cone and rod cells to^fettier. The rods respond to light intensity, rather 
25 than color. The mechanism of adaptation to cmihging intensity modeled by Dahari and 
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^Spitzeroperates in achromatic retinal ganglion cells 12. The RF of both chromatic and 
achromatic ganglieacells 12 includes both on-center and off-center receptors 14 of 
the center RF area, and receptoJ^J^iof the surround RF area. Hereafter, "center" and 
"surround" may be written as "cen" and "smcFS<§]Dectively. 

5 The responses of center area sub-regions 14 and surround area sub-regions 16 

are combined in ganglion cell 12, in one of two different ways, depending on the type 
of ganglion cell 12. As shown in FIG. 2, an "on-center" cell 20 responds positively 
("+") to signals of increased light intensity 24 (e.g. L + ) from the center area (or sub- 
region) of the RF, and negatively to signals of increased light intensity 26 (e.g. M") 

10 from the surround area (or sub-region) of the RF. An "off-center" cell 22 responds 
negatively ("-") to signals of increased light intensity 28 from the center area of the 
RF (e.g. L"), and positively to signals of increased light intensity 30 from the surround 
area of the RF (e.g. M + ). Cells 20 and 22 are also referred to as "opponent" cells. The 
physiological responses can be translated into mathematical expressions, and this will 

15 be done consistently throughout the following description. 

The "double-opponent" cells of the visual cortex (which have spatially larger 
RFs than the opponent cells), combine the responses of the on-center and off-center 
retinal ganglion (opponent) cells in a manner that resembles a mathematical spectral 
and spatial derivative. For example, FIG. 3 shows a double-opponent cell 40 that gets 
20 its own double-opponent center response Ld 0 -c ("double-opponent center" response, or 
"do-center" response) from a first group of on-center ganglion cells 42 (for example, 
L + M", the group that subtracts green surround responses from red center responses) 
whose RFs are located in a do-center area 44, and gets its own double-opponent 
surround response La 0 -s ("do-surround response") from a corresponding group of off- 
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center (e.g. L"M + ) cells 46, whose RFs are located in a do-surround area 48. The do- 
surround response is related to area 48. 



A "do-remote" area 50 has the shape of an annulus, concentric with that of the 
"do-center" and of the "do-surround" areas. Typically, the internal diameter of the 
5 "do-remote" area is equal to the external diameter of the "do-surround" and therefore 
does not overlap the "do-center" or the "do-surround"s. However, there may also be 
some overlap or some gaps between the do-remote and the do-surround areas 

One of the innovative aspects of the present invention is, that in contrast with 
US pat. 5,77 1,312, in the color contrast adaptation of the present invention the 
10 influence of a "remote" region affects the color contrast response of double-opponent 
(cortical) cells. 



cells) main groups of retinal ganglion cells 12 involved in color perception, 
corresponding to uVthree kinds of cone cells that respond to color analogously to the 

15 response of receptors (ro^sxells) 4 to intensity. An image is first processed in three 
most common color-coded champs in the retina (L + M~, M + L" and S + (L+M)~) to three 
activation-level maps of on-center o^onent P-LMS cells through different retinal 
layers originating in (receptors) cones 4 (FK^. 1). The off-center opponent cell types 
L~M + , M"L + and S'(L+M) + have similar RF strub^ure but with opposite signs. The 

20 input to the cones level is the spectral composition of the light reaching the retina, 
when illumination falls on a surface of the objects and is reflected from it. The field of 
view is mapped by the three types of cones, L, M and S. The quantum catches of the 
three cone types, L P j gme nt, Mpi gme nt and Spi gme nt, are typically expressed by an 
integration of the inner triple product of the spectral sensitivity of each pigment type 




Lere are six (or alternatively eight, including the yellow on- and off-center 
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with the reflectance properties and with the spectral and intensity absorption 
properties of the^^face, at any specific location in the image (Wyszeki & Styles). 
The responses of the three^JD^e-types, L CO ne, M con e and S CO ne, normalized separately to 
a range of 0-1, are typically expressfei^by a Naka-Rushton equation as a function of 
their inputs L p jg me nt, Mpi gmen t and S p i gme nt (Dalteri^and Spitzer; Wyszeki & Styles). The 
input red, green, and blue intensity values are treated a&vthough they are the responses 
of red, green and blue cone cells. 

The intensity values are transformed into a "response image" in a manner 
similar to that by which on-center and off-center retinal ganglion cell respond to 
10 inputs from cone cells and other cells in the retinal pathway. It is to be understood that 
in the present invention, in a manner analogous with that of US. Pat. 5,771,312, 
references to an "image" herein are references to values at pixels, or "pixel values", 
treated collectively as an array. Thus, the term "image" as used herein includes purely 
mathematical objects, and does not necessarily correspond to a physical image, 
1 5 although the original input images certainly do correspond to physical images. 

Th^s{n:st stage of the forward transformation is performed analogously to the 
tof> path of FIG/^Kof Dahari and Spitzer (notwistanding the fact that Dahari and 
Spitzer's paper does not contain color information), to provide, at each pixel, an 
opponent red center response (the^pixel value of a "response function"), an opponent 
20 green center response, and an opponeniblue center response. Optionally, a fourth, 
yellow opponent center response, which is calculated by simply adding the red and 
green responses, is obtained with the same forward transformation. In what follows, 
the fourth color (yellow) is optional in all transformations at th^ievel of opponent and 
double-opponent cells. A similar transformation is used to obtain reo^sgreen and blue 
25 surround responses. These surround responses are subtracted from rt^e center 
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r&spmises, in emulation of the action of "on-center" and "off-center" retinal ganglion 
cells, as folT^s^the green surround response is subtracted from the red center 
response; the red suiround re§j>€aise is subtracted from the green center response; and 
the red surround response and the greenSturound response are averaged to form a 
yellow surround response, which is subtracted fror?t*<the blue center response. 
Optionally, the blue surround response is subtracted from the yellowSfenter response. 

TheSgcond stage of the forward transformation deals with the simulation of 
double-opponent responses that emulate the action of (cortical) "double-opponent" 
color coded Pd 0 -LMSs cehs. and includes an additional remote adaptation. The color- 
coded double-opponent cells arcs^dapted (corrected) by a remote adaptation in a 
manner similar to a mechanism based c^psychophysical findings, as shown in Singer 
& D'Zmura, Vision Research, vol. 34, pp. 3H1-3126, 1994. Adaptation is also 
explained (but not in relation to color and do-remote arfe^) in Dahari and Spitzer. Note 
that "remote" adaptation refers to the effect of regions pen^heral to the "classical" 
RF. 

The third stage deals with the simulation of an "inverse" transformation", i.e. 
the transformation of the activity of these double-opponent "adapted" cells to a 
perceived image in a standard CIE notation (XYZ), in a RGB scale, in a Lu*v* scale, 
or in any other standard color scale. 

Regarding intensity contrast adjustment, the method includes similar forward 
and, optionally, inverse transformations, carried out by the emulation of the action of 
magno (M) cells. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The invention is herein described, by way of example only, with reference to the 
accompanying drawings, wherein: 
FIGTTls^^seheinatic cross section of a human retina; 

FIG. 2 is a schematic description of on-center and off-center opponent receptive 
fields; 

FIG. 3 is a schematic description of a double-opponent RF having a corresponding 

double-opponent remote sub-region; 
FIG. 4 is a schematic illustration of the operation performed to obtain a filtered 

opponent response. 

FIG. 5 is a schematic description of a double-opponent receptive field that includes the 
remote influences. 

FIG. 6 is an alternative schematic description of a double-opponent receptive field that 
includes the remote influences. 




1 5 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The mathematical simulation of the physiological processes related to color 
contrast perceptkm proceeds as follows: the spatial response profiles of the two sub- 
regions of each (on- N a*id off-) ganglion cell RF (P-RGC), the "center" and the 
"surround" regions, are preferat^ly expressed by a Difference of Gaussians (DOG). As 
in Dahari and Spitzer, the first step ofthe present invention is the transformation of each 
input image to an "output function" or "spech^d response function" G. For each color 
(red, green and blue), there is a center spectral response function G c (also written as L cen5 
M cen and S cen ) and a surround spectral response function G^(also written as L sr nd> M sm( j 
and S S rnd). Each of the three center response functions G c is computed by convolving the 
corresponding color image with a center local spatial Gaussian filter fcNjs shown in 
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ecju^tions 2 and 3 of Dahari and Spitzer, thereby producing a center smoothed image 
whose pikslvalues are the required G c 's. Similarly, each of the three surround response 
functions G s is com^ut^d by convolving the corresponding color image with an opposite 
sign surround local spatial Gattssjan filter f s? thereby producing a surround smoothed 
5 image whose pixel values are the requir^Gs's. Typical values of the radii p of the 
Gaussian filters are 0.5 pixels for the center filter aiiCkQ^-9 pixels for the surround filter. 

Optionally, the pixel values may be saturated beforeth^spatial filtering with the 
Gaussian filters: 

P := p/(p + Po) 

10 wherein p represents a pixel value, po is a saturation constant, and ":=" represents 
replacement. 

An emulated opponent response Ro P ( where R refers to L, M, S and optionally 
Y="yellow M ) is the subtraction between the center and the surround response values at 
a given pixel: for on-center cells the opponent response is expressed as L + M~, M + L" 
15 and S + (L+M)~, and for off-center cells as L"M + , M"L + and S"(L+M) + Opponent 
responses are calculated next. Specifically, for "on-center" opponent responses 

Rop(i+)= L + M = G Ci , - G s , m (1) 

Rop (m+)-M + L"= G c , m - G s j 

Rop (s+) =S + (L+M)= G c , s - (G s + G Sj ) ( , +m )/2 
20 Optionally, a yellow center response (Rop((m+i)+)) is also used. Similar subtractions are 
carried out to obtain "off-center" opponent responses L"M + , M"L + and S'(L+M) + (and 
optionally a yellow opponent response Ro P (( m +])-)) by subtracting the blue surround 
response from the yellow center response. 

The action of the RFs of double-opponent cells is emulated in one of two 
25 ways: 1) by assigning the emulated filtered response (explained below) (or response 
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in the specific case of feeding one opponent cell as the center of double opponent 
center receptive field region) of a first on-center opponent group (for example L*M") 
to a "do-center" response, and assigning the emulated filtered response of a 
corresponding, off-center opponent group (for example L"M + ) to a M do-surround M 
response, or 2) by using emulated filtered responses of "on-center" type opponent 
groups (or response in the specific case of feeding one opponent cell as the center of 
double opponent center receptive field region), for 

example assigning the emulated filtered response of an on-center L + M" group to a do- 
center response, and the emulated filtered response of an on-center M + L" group to a 
do-surround response. Hereinbelow, the "-" sign in all subscripts including two or 
more letters, for example in "do-c" of Ld 0 -c, means "of e.g., "L do . c ^ the red response 
of the double-opponent center". An exemplary mathematical formulation for the 
three-color double-opponent center responses is given by: 



center area 

M{do- c) = JJ/? * o^fi^ipc, yjdxdy 

center area 

S(do-c)= JJ/?* ol {s + )j '{^J&yfydxfy 



(2) 



center area 



a spatial and a spectrat^filtered response, similar to Ro P (i+) 5 Rop (m+), Ro P (s+) (and 
optionally the fourth "ydlow^N^op^i+m^) but related to different areas than the 
corresponding Ro P s. Ro P and R* op and similar pairs below are also referred to herein 
as, respectively, "non-filtered" and "filtered" responses (in spite of the fact that the 
"non-filtered" values are obtained of course with opei^tions involving filters). The 
relationship between non-filtered and filtered opponent responses is illustrated 
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schematically in FIG. 4, where an examplary (non-filtered) on-center opponent 
response Ro P (xi,yi) is calculated for pixel 62 centered at (xi,yi) using eq. 1, i.e. 
subtracting the value at (xi,yi) of a surround response function G s (xi 5 yi) from a center 
response function G c (xi,yi). As explained above, both G s (xi,yi) and G c (xi,yj) are 
5 functions obtained by convolutions involving local filters. A filtered response 
Ro P *(xi ? y]) to be used in eq. 2 for the computation of the do-center responses for the 
do-cell centered at (xi,yi), (and thus associated with Rop(xi,yO) is calculated in 
general by applying the local filters in the convolutions yielding G c and G s at different 
(than xi,yi) locations 62. The "filtering" operation for obtaining filtered responses 

10 thus involves moving the filters over the relevant area of the integration and obtaining 
for each pixel a function value. The value of R* op is in general different from that of 
Rop 5 except for the special case when both are calculated at the same (x,y). 

Optionally, a fourth yellow do-center response is obtained as ((M+L)-S)d+ C en 
by using R* op ((i+m)+. Thus, the responses L(do-c), M( do - C ) and S(do-c) (and optionally 

15 (M+L)-S)a+cen) are the products of convolutions of corresponding filtered opponent 
responses with a center spatial weight function f( C -do)- f(c-do) is preferably a radially 
symmetric Gaussian exponentially decaying function, with a radius of 1-6 pixels. 

The do-surround responses are computed in a similar fashion. An exemplary 
mathematical formulation for the three color double-opponent surround responses 

20 (and optionally a fourth, yellow (L+M) ( do - S ) ) is given by: 



ffi*of(Hf« s ^x,y)(briy 



(3) 




surroundirea 




surroundirea 
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where R* op (i-) 5 Ro P *(m-)» and Ro P *( S -) ( and optionally a fourth "yellow" R* op ((i+m)-) are 
off-center filtered responses similar to Ro P Ro P (m-)> Rop(s-) (and optionally Rop((i+m)-) 
but related to different areas, as explained above. These off-center filtered opponent 
responses are convolved with a surround spatial weight function f( S -do> f(s-do) is 
5 preferably an radially symmetric Gaussian exponentially decaying function with a 
radius of 9-31 pixels. 

As mentioned above, a key innovative element of the present invention is the 
addition of a new, "do-remote" signal needed for the adaptation mechanism of color 
contrast, for each double-opponent cell. This "do-remote" signal represents the 

10 peripheral area that extends far beyond the borders of the classical RF of the Pd 0 -LMS 
double-opponent cell. The four do-remote signals (L( do -remote) 5 M( do -remote) 5 S (do-remote) 
and (L+M)( do -remote) that emulate the effect of the remote areas of the double-opponent 
cells are defined in eq. 4 as the convolution of a corresponding filtered double- 
opponent center response with a remote spatial weight function (representing spectral 

15 content) f r applied to do-remote area 50, for each pixel in the remote area. 

L(do - remote) = * (Jo _ c) f r (x, y)dxcfy (4) 

remote-area 

M(do - remote) = JJm * {do _ c) f r (x, y) dxdy 

remote-area 

S(do - remote) = ||s * {d0 _ c) f r (x, y)dxdy 

remote-area 

Y(do - remote) = (L + M)(do - remote) = || (L+M) {do c) ^ ^ y)dxdy 



remote-area 



where L* (do -c), M* (do -c), S* (do -c), and (L+M)*( do -c) = Y* ( do -c) are double-opponent 
center filtered responses acting on the do-remote areas, related to the responses of eq. 
2 in the same fashion as explained above for opponent responses. In other words, the 
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same calculations that were carried out for opponent RFs are done here, but in relation 
to each pixel that belongs to the do-remote area. 

f r is preferably an exponentially decaying spatial-weight function, an example 
of which is shown in Eq. 5. Alternatively,/, can be a Gaussian function, or any other 
decaying function. 



exp 



(5) 



Krembte is a constant which defines the slope of the weight function, and A remo te is a 
factor of normalization to a unit: 

^ (6) 



JJ exp 



remote-area 



The next step is obtaining the responses of the double-opponent RFs that are fed by 
10 both the do-center responses (eq. 2) and by the do-surround responses (eq. 3). The 
double-opponent responses or "output functions", before the adaptation stage are 
given by: 

Ldo = L(do - c) — M {do - s) (7) 
M do = M {do - c) — L (do - s) 
Sdo = S (do - c) — (jL + A/ )(c/o - a-)) 

Optionally, a fourth (yellow) double-opponent output function is obtained by 
subtracting the blue do-surround response from the yellow do-center response. The 
15 do-remote signals of eq. 4 can alternatively be built from filtered double-opponent 
responses L*(d 0 ) 5 M*( d0 ), S*(do) and Y* ( do) = (L+M)* (d0 ) (related to the do-responses in 
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eq. 7 as explained for the opponent case above) through appropriate convolutions with 
the remote spatial weight function, instead of using the double-opponent center 
-^jUtered responses L*(do-c) 5 M* ( do-c) 5 S*( do -c) and Y* (do -c). 

^2^\ J ^^an be understood from the subscripts "c" (center), M s M (surround), "do-c" 
5 (double-oppohent center), "do-s" (double-opponent surround) and "do-r" or "do-remote" 
(double-opponenwemote), the various receptive fields ("RF") and spatial convolutions 
extend over differenr^iumbers of pixels. FIGS. 5 and 6 show two different schematic 
descriptions of double-opponent RFs which include the remote influences. FIG. 5 shows 
an on-center cell 100 with a\enter RF 102 (pixel A) and a surround RF 104 (pixels B). 

10 The on-center cell is suiroundeckby 9 off-center cells 106, each with a center RF 108 
(pixel C) and a surround RF 110 (pixels D). Cell 100 and cells 106 are "opponent cells", 
and together constitute a double-opponetat ("do") cell 112, on-center cell 100 constituting 
a "do-center" RF, and off-center cells 106\x)nstituting a "do-surround" RF. Do-cell 112 
has a do-remote RF 114, which includes on-oenter opponent cell RFs 116, each with a 

15 center RF 118 (pixel E) and a surround RF 120 \pixels F). For simplification purposes, 
only nine such RF 116 are shown in the upper rignt corner of do-remote RF 114, with 
the understanding that they are repeated in the res^)f the (unmarked) pixels. The 
different marking of pixels in cell 100 and RFs 116 Indicates the different values 
resulting from the application of different spatial weight functions at those locations. In 

20 the example of FIG. 5, the different receptive fields extend overdifferent areas: for both 
on-center and off-center opponent cells, the center RF includes None pixel, while the 
surround RF extends over the 8 pixels immediately surrounding each ctenter. For double- 
opponent cell 112, do-center RF 100 extends over 9 pixels and the oo-surround RF 
(which includes the nine cells 106) extends over the adjacent 72 pixels. DVremote RF 

25 114 extends over the 648 pixels surrounding do-cell 112. It should be cleak that the 
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cKff^rent areas over which the spatial weight functions are applied in the convolutions 
can also partially^vSriap^orbe separated by gaps. For example, cells 106 whose RFs 
compose the do-remote RF and who areu§5d-itLComputing the do-remote signal can 
partially overlap, or can be separated by gaps. 
5 FIG. 6 shows an alternative double-opponent RF description. The main 

difference vs. FIG. 5 lies in the composition of the do-remote RF: instead of (as in FIG. 
5) being composed entirely of on-center cell RFs, a do-remote RF 130 is composed here 
of a group of do-cell RFs 132, each RF 132 composed of a center "on-center" opponent 
RF 140 (pixels E and F) and nine surrounding "off-center" opponent RFs 142 (pixels G 
10 and H). For simplification purposes, only one of the nine RF 132 is shown in detail. The 
pixels in the on-center and off-center RFs of the do-center, surround and remote regions 
are again marked differently to indicate the effect of the convolutions with different 
spatial weight functions. The do-remote area extends over the same number of pixels as 
in FIG. 5 

15 The descriptions in FIGS. 5 and 6 are for illustration purposes, and the areas 

given are typical. However, it is to be understood that each area mentioned can extend 
over a range around the typical value cited. Thus, an opponent surround RF (e.g. 104) 
may include the 16 pixels immediately surrounding the eight pixels mentioned above, 
thereby extending over a total of 24 pixels. The RF of a double-opponent center typically 

20 extends over about 3 by 3 pixels, see FIG. 6. The do-surround RF may extend over an 
area defined by another 11 pixels extending each way in both x and y (i.e. bound 
externally by a boundary of 31x31 pixels) while the do-remote RF may extend an 
additional 10 pixels in each direction (i.e. bound externally by a boundary of 51x51 
pixels, and internally by the 9x9 to 31x31 boundary). The domains of the convolutions 
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for computing the different terms roughly overlaps the corresponding RF areas. At the 
boundaries of the images, all convolutions use periodic boundary conditions. 



Adaptation 

5 The adaptation for color contrast correction is carried out on the double-opponent 
responses of eq. 7, in a manner analogous to the adaptation procedure in Dahari and 
Spitzer, by using the do-remote signals. In the form of the present invention that is 
applied to digital video photography, time variation is taken into account. The output 
functions G (where G stands for Ld 0? M do and Sdo> the double-opponent output 

10 functions in eq. 7) are now functions of time t 5 specifically the sampling time of the 
images. In addition to the spatial and spectral weight functions which are applied 
above to entire images at a single value of t ? temporal filters, as defined below, are 
applied pixel-wise to these functions of t. Two kinds of temporal filters are used: 
temporal, preferably low pass filters, which are independent of the functions to which 

15 they are applied, and adaptive temporal filters, which depend, at any given time t, on the 
prior history of the functions to which they are applied. 

For each color, the present invention computes an adapted (or "corrected") 
double-opponent response R<do-a)(G,t) at each pixel by using a Naka-Rushton-like 
equation: 



20 R max is a maximum response that is used as a normalization factor. For the remote 
adaptation of each double-opponent RF region, an adaptive function Gb(t) analogous to 
Dahari and Spitzer's adapting component Gb(t), is used in the computation of the 
semisaturation term a,: 




(8) 
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a(t) = aGb(t) + P 



(9) 



In the application of the method to still photography (no time dependence) i.e. at steady 
state (at t=oo), function Gb is time-independent and is given (using the red double- 
opponent response as an example) by: 



5 Similar equations are used for the other colors response components (green, blue and 
optionally yellow). Unlike Dahari and Spitzer, the Naka-Rushton equation of the 
present invention is applied to the double-opponent RFs (emulating the cortical color- 
coded cells), and the adaptation factor (the "remote influence") is applied to the semi- 
saturation term in the denominator of eq. 8. More generally, for example in video 

10 applications, the adaptive function is computed by convolving the corresponding output 
function with corresponding adaptive temporal filters as shown in equations 7 and 8 of 
Dahari and Spitzer, and in the boxes labeled "dynamical adaptive filter" of FIG. 1 of 
Dahari and Spitzer. The temporal filters are applied pixel-wise and locally to the entire 
image. What makes these filters adaptive is that the associated time "constants" actually 

15 are functions of both time and the prior history of G. The most preferred form of the 
function ib that describes the decay of the adaptive filter is: 



In this expression, x m is the maximum expected value of Xb(t), G(t) is the output function, 
after convolution with the corresponding low-pass temporal filter; Gb(t) is the adaptive 
20 function, i.e., the output of the convolution, at times prior to the time at which the 
convolution presently is being computed; and G n is a normalization constant. Because 
the adaptive filter is causal, it is well-defined despite being defined in terms of its own 
output. 



Gb(t^co) — L(d 0 -remote) 



(10) 



Xb(t) 



= x m /( 1 + abs( G(t) - G b (t) )/G n ) 



(11) 
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The change in a gain control effect is achieved by the curve-shifting 
mechanism of the "response vs. log chromatic (or achromatic) contrast illumination" 
curve, similar to the "response vs. illumination" curve shown schematically in FIG. 6 
in Shapley & Enroth-Cugell, Progress in Retinal Research, vol. 3, pp. 263-346, 1984. 
5 The gain control effect of the present invention can be seen as a "second order" effect, 
in contrast with the first order effect embodied in FIG. 6 of Shapley and Enroth- 
Cugell. The adaptation is reflected in a shift of the response curve as a function of 
time. Consequently, each time a new range of input intensities to a color channel is 
viewed, given that the change from the previous stimulation is sufficiently large, the 
10 curve shifts, bringing the system to a new adaptation state. This curve shift takes place 
according to the temporal filter, causing an apparent decaying function of the 
response. 

The significance of applying a variable dynamic temporal component, aside 
from corresponding to known physiological phenomena (e.g. as explained in Dahari 
15 and Spitzer and the references therein) becomes evident when observing the 
dynamical color contrast adjustment algorithm performance in the absence of such a 
component. In such a dynamical case, the rate of adaptation changes as a function of 
signal magnitude of each contrast channel. 

One immediate result of the dynamical properties of adaptation mechanisms in 
20 video applications, (as also explained in U.S Pat. No. 5,771,312) is that when a 
constant image is viewed over a certain time period, its color contrast appears to 
decay with time. 

"Inverse transformation" or transformation of the Pdo-a LMSs activity to a 
perceived image 
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In \rder to perform contrast on real images it is necessary to inversely 
transform the Adapted (corrected) double-opponent responses (at any location in the 
image) into a perc^ved color. Several different inverse functions, based on different 
assumptions, can be uted for this purpose in preferred embodiments of the present 
5 invention. The calculatecKperceived color contrast is the color contrast that would 
stimulate the triplet of Pdo-a \MS cells to the same responses, with a uniform non- 
contrast surface present in thei\ "remote" areas. The rationale that has lead to this 
definition is the assumption thatNthe visual system interprets color contrast of an 
object in a way similar to the one described above. The following equations describe 
10 the emulation of the three main double-exponent and opponent color-coded cells, but 
alternatively all these equations can be applied also to the yellow double-opponent 
and opponent RFs. In general, the steps in \he inverse transformation preferably 
include: using the adapted double-opponent responses to calculate a new value of a 
double-opponent response ("new do-response"), assigning the new value to the 
15 double-opponent center ("new do-center response") thiife obtaining a new value of an 
opponent response ("new opponent response"), and assigning this new opponent 
response value to get a new opponent center response, thus\eturning to eq. 1, which 
now has "new center" values. 

A key assumption in performing an inverse function as described above is that 
20 the contrast in the "do-remote" area is equal to that in the double-opponent area. This 
is expressed as 



(12) 



j j -is \yf o 

do-a ~ do -remote do-a do -remote do-a do -remote 
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Inserting for each color the do-remote expressions of eq. 4, and rewriting eq. 8 
for each color channel under assumption (12) yields: 



^do-a 



L do (13) 



4*. + \\ L *jo-c f r <<x,y)dxdy 



remote 

M, - 



1 do- 



° M do + jjM* do _ c f r (x,y)dxdy 



remote 

s _ S 

^do-a — 



S*+ j\s* do _ c f r (x,y)dxdy 



By assuming for each color channel that the semisaturation term expressed by 
5 the integral in the denominator equals approximately the adapted double-opponent 
response (e.g. as an alternative for the red Ld 0 -a~JJL*do-cfr(x,y)dxdy) we get 



T t_ i^do-a ) 



(14) 



(s Y 

o t V do-a ) 

1 ^do-a 

where Ld 0 ' 5 Mdo' and Sdo' (and optionally Y do) are the new double-opponent responses. 
Ldo' ? Mdo* and Sdo' (and optionally Y'd 0 ) are now applied to the corresponding double- 
opponent center color channel, i.e., L'do — L'd 0 -c> M'd 0 = M'd 0 -c, S'do = S'do-c, (and 
10 optionally Y' d o = Y' do -cX at each pixel, creating for each color a difference ("change") 
between the original double-opponent center response (Ldo-o Mdo-c, Sdo-c ? eq. 2) and 
the new do-center response (L'do-o M'do-c, S do-c) 
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U da _ e (x,y) = L db _ e '{x,y)-L dB _ c ( X ,y) = - L *°' L - do ° 5) 



bMjo-c (*, y) = M do _ c ' (x, y) - S Jo . c (x, y) = M "° ' Md ° 



*S Jo -c (*, y) = S Jo _ c ' {x, y) - S do _ c (x, y) = 



fc- d oix,y)\\dxdy 

cen 

fc-doi^y)\\ dxd y 



Remembering that do-center responses are equivalent to opponent responses 
(connected by a convolution with a weight function) the new do-center response 
values are equivalent to new opponent response values, in the special case when the 
do-center RF equals the RF of one opponent cell. 

In order to find the new opponent center responses, i.e. the new responses of 
the cones, the AL(do-c), AM(do-c), and AS(do-c) differences are applied to the 
opponent center, preferably using the following expressions: 

r cen < x > = <* AL (do-c) (*> + bL cen (*> >0 + cM (*> srnd 
M% cen (*>^ = aAM (do-c)( x > + bM cen( x > + cL( < x > yhrnd 



S 'cen (*> y) = aAS (do-c) (*> y) + bS cen <*> ^) + C ( L + M * x > * W 

where a ,b and c are constant parameters that can be identical for all color channels, or 
alternatively are different for each color of the center sub-region. Typical ranges for a, 
b, and c are a=l-2, b= 0-1 and c=0-2. The outputs of eq. 16, Lewi, M' C en, and S' ce n are 
the new values of L 5 M and S. Alternatively, the new cone response values can be 
obtained with a multiplication operation instead of a summation operation, for 
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example using the following type of expression (for the red channel, with similar 
expressions for the other color channels): L' cen (x,y)= aL ce n(bL'do-c(x 5 y) + cL cen (x,y) 
+dM srnd + e) where a, b, c, and 8 are constant parameters with typical value ranges of: 
a=l-5; b=l-5; c=0-2, d=0-5, and 8, a constant which is normally 0. For specific uses 
such as correcting JPEG images, "a" in eq. 16 and in the multiplication operation 
above can be assigned negative values (e.g. -1 to -5). Thus, modified new cone 
responses can be obtained either by an additive operation, i.e. by adding a change to 
the original value, or by a multiplication operation as shown above. Another 
alternative for obtaining the new opponent center (i.e.cone) responses through an 
additive operation as in equation 16 is (for the red channel as an example) 
L'cen(x,y)= bLdo-c(x,y) + cL cen (x,y) +dM srn d + 8, where the typical value ranges are: 
b=l-30, c=0-l, d=0-5, and s= normally zero. For specific uses such as correcting 
JPEG images, "b" in the additive operation above can be assigned negative values 
(e.g. -1 to -30). 

Finally, a simple inverse transformation from the LMS scale to the CIE-RGB 
(either directly or by going through the CIE-XYZ scale) can be performed. It is 
important to note that the method of the present invention is independent of the color 
scale used, and that there are other scales, such as the non-linear Lab scale, to which 
the method of the present invention is equally applicable. 



Intensity Adaptation 

Intensit^adaptation or contrast adjustment on the achromatic information in a 
coior image is connectecP&>tfcie color contrast adaptation described above, in the sense 
that it uses a similar physiological baSk^for the algorithm and a similar sequence of 
steps, except that in one alternative case ("case r*)4he adaptation is performed on the 
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center and surround areas of the opponent RFs, and not on the double-opponent RFs as 
in the cbk>r contrast case. Another difference vs the color contrast adaptation is that in 
the intensity ostntrast adaptation of case 1, an inverse transformation may not be applied. 
Intensity adaptation can be applied independently on achromatic images (non-color or 
5 black-and-white pictures) oiN^n the intensity domain of color images. For intensity 
adaptation, the method emulates twcNjmes of ganglion cells of the magno pathway, the 
on-center and off M-opponent cells. In a\preferred embodiment of the intensity 
adaptation, the remote adaptation is applied to the ceht^r sub-region of opponent RFs (or 
alternatively to the center and surround sub-regions of the op^x^nent RFs), unlike in color 

10 contrast, in which remote adaptation is applied to double-opponenNRFs. Alternatively, 
achromatic contrast adjustment can be performed on double-opponent RFsHpo. 

The information contained in a color or black-and-white scene is processed 
through rod cells 10 (FIG. 1) or through the three cone types. In the first stage in the 
"achromatic contrast" or intensity adaptation method (algorithm), an emulation of this 

15 action is done under the assumption that the three cone types absorb the same amount 
of light. Alternatively, the starting point of the forward transformation in this case, as 
in the color contrast case, can be an image acquired using any of the means described 
above. Regardless of the color scale in which the image is obtained - RGB, XYZ, 
(Yxy), Lab, or any other CIE scale - the intensity value at each pixel can be obtained 

20 by known mathematical manipulations. Thus, the starting point is a known intensity 
value at each pixel. 

The steps in the intensity adjustment method start in a similar fashion with that 
in the color contrast adjustment method, with the distinctive difference that one 
emulates rod cells and not cones. Thus, if the input is the intensity spectrum of a 
25 scene, an emulated pigment response Ypj gme nt is given by 
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visual range 



where I ro d(^) is the spectral absorption function of the rods as a function of the 
wavelength X; 9?(A,) is the reflectance property of the surface at any specific location 
in the image; and is the spectral composition of the illumination falling on this 
location. 

5 In a similar fashion, the emulated responses of the rod- types at each pixel (or 

the intensity value at each pixel) R ro d are expressed by a Naka-Rushton equation as a 
function of their inputs Y p i gme nt ? 

Y» (18) 

D — pigment 
rod yn n 

1 pigment °(AW) 

Here too, in the case that the intensity information is extracted from rods, parameters 
n and ctnr are similar to those related to cones (see eq. 2). Alternatively, eqs. 17 and 
10 18 are not relevant and are not used if the starting point is an existing image. 

Continuing along the lines of the color contrast equations, a "center" response 
before the adaptation is now defined as 

center area 

whe$£ the RF region is^pcgferably circular, f c is preferably a Gaussian decaying 
spatial-weight function, analogous to^thotin eq. 4. The "surround" response, which 
15 now represents the surround sub-region of the RFcrfs^n M-RGC, is defined as 
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As mentioned already, a key difference in intensity contrast vs. the color 



5 contrast adaptation is that in the present invention a "remote" component needed for 
the intensity adaptation is computed already at the opponent RF stage for both 
"center" and "surround" regions of the on-and off-center RFs. The "remote" signal (or 
"response") representing an adaptation component derived from peripheral regions, 
i.e., areas that extend outside and beyond the classical RF ("remote adaptation"), is 
10 calculated separately for each of the sub-regions (or alternatively applied only to the 
"center" area of the on-center and off-center cell types) using an equation of the type: 

G_ = \\G * cer,f r {x,y)dxdy (21) 

remote area 

where G* ce n is a filtered response (with a definition similar to that in the color contrast 
case) similar to the center response in eq. 19, but related here to the remote area 
through the convolution with a remote spatial weight function f r . Alternatively, an 
15 additional remote signal can be built with a surround filtered response G* srn d similar 
to the surround response of eq. (20), but related here to the remote area (by 
convolution with a remote spatial weight function f r , and using an integration similar 
to that in eq. 21). 

The mechanism for the intensity adaptation application can now be calculated 
20 in one of two ways. One way (Case 1) is to apply the adaptation and the remote area 
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influence at the emulation stage of ganglion RF regions (eqs. 20, 21) or alternatively 
(Case 2) at the stage of the cortical level, in which case the adaptation is applied on 
the double-opponent RFs (similar to the color coded cells case). 

Case 1: The M-RGCs 'center' and 'surround' sub-regions are adapted 
separately (the values of the parameters are based on electrophysiological findings, in 
a similar manner to that found and analyzed for M-RGCs by Shapley & Enroth-Cugell 
and Dahari & Spitzer). As described in these references, and as for color contrast, the 
response of an "on" and "off" M-RGC to a stable visual stimulus at t=oo (the "static 
case") is: 

(22) 



Rop (/ = oo) = 



smd 

G cen + o" cen G smd + <y srnd 



10 Where a is an adaptation factor which depends on the stimulation (here, at t=oo) of 
both the sub-regions and that of the remote area. In other words, the adapted opponent 
response is obtained by a subtraction of the adapted surround response from the 
adapted center response. The adaptation factor of each sub-region has separate local' 
and 'remote' components and consists of the same spectral properties as the sub- 

15 region. For example, for the "center" sub-region, the adaptation factor is given by: 

^ cen -cr cenJocal+ a c „ e (23) 
where cr cenJocal = a cen G cen + p cen 

QYlCi & cen, remote ^cen remote 

The adaptation factors for surround sub-regions are similarly defined, and include 
both local and remote adaptation components, a, p and c for both center and surround 
adaptation factors are constants with different values for the 'center' and for the 
'surround', a and J3 are determined as in Dahari and Spitzer, while c in eq. 23, the 
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weight function of the remote adaptation component, typically ranges between 1-30 
for both center and surround adaptation factors. The method can be alternatively 
applied without the surround signal and its remote mechanisms, i.e. using only the 
first term in eq. 22. Furthermore, the method can be alternatively applied without 
considering the local adaptation, i.e. using only the second term in eq. 23. 

If one defines a= a +1 and b=|3 and substitutes the terms of eq. 23 in eq. 22, 
the emulated responses of the on-center and off-center cells M-RGC are written as: 

G G , (24) 

D cen srnd ^ / 

K op — 

^ cen G cen ^ ^ cen ^ cen ^ remote ^ srnd ^ srnd ^ srnd ^ srnd ^ remote 



For the dynamical case, one can follow a similar sequence to that expressed 

10 by eqs. 9-1 1, to be used here however for achromatic intensity. 

. Case 2: In the alternative embodiment of performing the adaptation for 
intensity contrast by applying it at the double-opponent stage, eqs. 22-24 are not 
implemented. In general, the steps of the method for intensity adaptation at the 
double-opponent stage follow those of the color contrast adaptation. One difference 

15 vs. the color contrast case is that the input of the do-surround sub-region is optional, 
and intensity adaptation can be performed without it. In analogy with eq. 1, the 
opponent cell's response, Ro P is the subtraction between the center and the surround 
output function values of on-center and off-center cell types: 

a, - c c (25) 



32 




where G cen is equivalent to G c and G srn d is equivalent to G s , both however being 
applied here for the cells that code the intensity domain. The off- center cells have 
similar RF structure but with opposite signs. 

In analogy with the color contrast case, the achromatic double-opponent center 
5 signal (do-center response) is built from the filtered responses of the opponent on- and 
off-center cells: 

center area 



where R* 0 p is a filtered response similar to Ro P but related to a different area (as 

explained in the color contrast case), an area defined in the integration. The typical 
size of the do-center contains 1-16 pixels in the simulations applications. 

10 The' do-surround' signal is similarly defined as 

Qdo-s) = ffi*of(do- S )(x,y)cbdy (27) 

surrowwl_area 

The do-remote signal, in analogy with eq. 4 is given by 

GU™*= ^* {d ^f r {x,y)d^y (28) 

remote area 



The double-opponent responses (the responses of the on-center and off-center 
double-opponent M-RGCs), before the adaptation stage, are given by: 



33 



34 



Rdo 



— G (do - c) — G (do - s) 



(29) 



The do-remote signal can alternatively be built from the filtered response R*do of the 
double-opponent RF. The adaptation, for still and video images, as well as the inverse 
transformation are then performed analogously to those in the color contrast case, 
following the same sequence of steps, except for the fact that it refers here to 
achromatic information, and that, as mentioned above, the intensity adaptation can be 
performed without the input of the do-surround response. 

All publications, patents and patent applications mentioned in this 
specification are herein incorporated in their entirety by reference into the 
specification, to the same extent as if each individual publication, patent or patent 
application was specifically and individually indicated to be incorporated herein by 
reference. In addition, citation or identification of any reference in this application 
shall not be construed as an admission that such reference is available as prior art to 
the present invention. 

While the invention has been described with respect to a limited number of 
embodiments, it will be appreciated that many variations, modifications and other 
applications of the invention may be made. 
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